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ABSTRACT 


The frequency response analysis of several constrained 
viscoelastically damped plates were performed with 
modifications to account for the frequency dependent 
properties of viscoelastic materials (loss factors and shear 
modulus) using MSC/NASTRAN. A modal strain energy 
approximation method was also applied to each plate. 
Comparisons of the frequency response analysis to the modal 
strain energy method were made for each plate structure. 
The experiments were performed on the plates to establish 
the base line data for the comparisons. 

The modal strain energy method frequency response plots 
showed good agreement with the frequency response analysis 
method plots for all plate structures analyzed. Plate 
structures with fewer natural frequencies in the frequency 
range investigated showed better agreement between the two 
prediction methods than the structures with more natural 
frequencies. The experimental data for the plate analyzed 
indicate less damping than calculated by either of the 
finite element prediction methods. The experimentally 
determined natural frequency modes were well approximated by 


both of the finite element prediction methods. 
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Lee ODUC TION 


Advances in the material sciences area have resulted in 
metals with high strength to weight ratios. The application 
of such material can result in structures with low levels of 
damping. Exposure of these structures to high levels of 
excitation creates a need for damping techniques which are 
both inexpensive and effective. Introduction of a 
distributed viscoelastic material with an incorporated 
constraining layer has gained widespread acceptance as a 
means of reducing vibration response. Many finite element 
codes have been developed which incorporate frequency 
prediction methods for use in design applications, reducing 
the need for experimental analysis. However, accurate 
behavioral predictions of applied viscoelastic material have 
remained costly, due to the frequency and temperature 
dependence in viscoelastic materials. 

The response of the simplest structures can be 
represented by continuous classical methods but the finite 
element approach remains the method of choice when complex 
physical systems must be analyzed. The use of a finite 
element method to find the system response requires 
modifications to capture the environmental dependencies of 
the applied viscoelastic layer. One standard method of 


conducting a frequency analysis of dampened structures is 


through repeated application of the frequency response 
analysis over incremented steps in frequency. A second 
approach, using modal strain energy approximations in 
conjunction with the finite element method, has proved 
effective in predicting excitation responses of structures 
with a constrained viscoelastic layer. This paper 
investigates the two prediction methods described above as 
applied to plate structures. MSC/NASTRAN was chosen for the 
analyses because of its widespread availability and previous 
use in this area. The prediction methods were applied to 
two plate structures, one with the damping material having 
edge exposure to the environment and a similar plate with 
the damping material protected from the environment and 
possible corrosion, through the use of an edge weld. The 
edge-welded plate was evaluated before and after welding to 
determine possible damage to the viscoelastic layer. The 
prediction methods were evaluated for ease of application 
and prediction accuracy as compared to the experimental 


data. 


II. THEORY AND BACKGROUND 


Continuous structures in the physical world possess 
distributed characteristics of mass, damping and stiffness. 
(Ref. lspp. 140-141] Structures with such characteristics 
are said to contain an infinite number of degrees of 
freedom. The classical vibration analysis of such a 
continuous system involves the construction of a discrete 
mathematical model with a finite number of parameters to 
approximate the overall system response. Construction of a 
model of N discrete elements allows the discretized 
continuous system to be represented by the coupled matrix 


equations: 


[M]{q(t)} + [C]}{q(t)}) + [K](aq(t)}) = (Q(t)} (2) 
Linear transformation and normalization of the undamped 
portion of Equation (2.1) and assumption of damping 


approximated by an appropriate diagonal matrix gives the 


uncoupled equations: 
Ca(t)) + C 2cwn 1€4(t)) + Log JC a(t)) = (1(t)) (ene) 


(q(t)} = [o]{ a(t) } (2.3) 


where: 


displacement vector 


{q(t) } 


[o] = mode shape matrix 

[ we ] = vector of natural frequencies 
n 

C = damping factor 


fey } transformed matrix of applied loads. 


Application of modal analysis through developed finite 
element codes have resulted in readily obtainable solutions 
to Equations (2.2) and (2.3) for complicated, lightly damped 
structures. [Ref. 2:Secs. 203, 4.3527 Ref. cor 2] 
Recently, viscoelastic materials have been added to 
structures to improve their overall damping characteristics. 
However, the addition of these viscoelastic materials to 
design structures introduce both frequency and temperature 
dependent variables that must be incorporated into the 
finite element analysis. [Ref. 4:p. 2] 

Viscoelastic damping has been observed in many polymeric 
materials and has been directly attributed to the relaxation 
and recovery of long carbon atom chains of the polymer after 
deforming. [Refs “> Joo? | Because of their effect on 
molecular motion, temperature and frequency influence these 
"relaxation and recovery" properties. Thus, a complete 
understanding of the temperature and frequency dependent 


behavior of the material is essential. 


Moderate changes in temperature or large changes in 
frequency result in definite predictable alterations in the 
viscoelastic material properties of shear modulus of 
elasticity and loss factor. Actual quantitative changes in 
the damping properties vary substantially with the chosen 
material. Qualitative changes are categorized into three 


descriptive regions; the glassy regions, transition region, 


and the rubbery region as shown in Figure 2.1. Ina given 
viscoelastic material, all three regions are typically 
contained in two to three decades of frequency. [{Ref. 5:p. 
91] 


Damping at relatively low temperatures for the given 
material is in the glassy region and exhibits a low loss 
factor and high shear modulus of elasticity. As the 
material temperature iS increased the shear modulus begins 
to drop slowly as the loss factor rises. Surf rerent 
temperature increases allow the damping material to enter 
the transition region, characterized by a rapidly decreasing 
shear modulus value. The loss factor reaches a peak value 
at the maximum rate of change in shear modulus value. 
Further increases in temperature result in the material 
entering the rubbery region with a corresponding decrease in 
loss factor and shear modulus values. While the above 
discussion is descriptive of varying properties with 
temperature increases, it also applies to  £frequency 


decreases. Examination of several materials shows a 25 


degree fahrenheit increase being roughly equivalent to a one 
decade frequency decrease. For example, a 25° F increased 
temperature is roughly equivalent to a 10 kHz to 1 kHz ora 
1000 Hz-100 Hz drop in frequency. [Ref. 6:pp. A1-A21] 

The high sensitivity of viscoelastic materials causes 
the matrix values of Equation (2.2) to change with only mild 
variations in excitation and environmental conditions. This 
results in finite element solutions that are only applicable 
over a narrow range of frequency or temperature. A solution 
to Equations (2.2) and (2.3) which reflects varying matrix 
values may thus be obtained by the segmentation of the 
equations into incremental modal frequency response finite 
element analyses steps. This composite approach gives good 
results but in practical applications may become 
prohibitively expensive and time consuming, especially when 
several design alternatives are under consideration. 

An application of the modal strain energy method (MSE) 
by Johnson and Kienholz [Ref. 7:pp. 5-12] is suggested as an 
alternative solution to Equations (2.2) and (2.3) by 
imposing frequency dependency on the frequency response 
through post analysis data manipulation. In the MSE 
approach, the MSC/NASTRAN code is used to construct a finite 
element model with constant, midrange values of loss factor 
and shear modulus of elasticity. NASTRAN modal analysis is 
applied and the mode shapes and natural frequencies are 


extracted and inserted into Equations (2.2) and (2.3). The 


fraction of elastic strain energy attributed to the 
viscoelastic core at each natural frequency is also obtained 
from the modal analysis and the damping material loss 
factors interpolated at each natural frequency. The matrix 
[2ctw,] of Equation (2.2) is then replaced with int are 
when Equations (2.2) and (2.3) are modified CRN the 
interpolated loss factors, the frequency response solution 
becomes frequency dependent. The MSE method thus gives 
accurate solutions over a large frequency range while 


remaining cost effective. 


G, Or Loss Factor, n 


Shear Modulus, 


Constant 
G Frequency 


| Glassy Region Transition Region Rubbery Region 





Tenpemcre Wise 


Figure 2.1 Variation of Shear Modulus and Loss 
Factor with Temperature [Ref. 5:p. 90] 


III. MATERIAL AND METHODS 


A. MODEL DESCRIPTION 

Three different plate structures were analyzed during 
this study in order to provide a basis for the comparison 
with the finite element results. These plate structures 
consisted of a simple plate, a milled plate with an 
unrestricted constraining layer and finally a milled plate 
with a restricted constraining layer. The simple plate was 
a 14 x 30 inch aluminum rectangular plate, 3/8 inch thick 
(Figure 3.1). Dimensions were chosen to give approximately 
15 resonant frequencies under 1000 Hz. The milled base 
plate and its dimensions are shown in Figure 3.2. All 
plates were machined from standard 6061 T6 aluminun. The 
use of a damped plate structure in a corrosive environment 
requires design modifications to protect the viscoelastic 
material. The welded plate structure of Figure 3.2 is such 
a design using the same material and overall thickness as 
the simple plate. This milled plate was subjected to 
vibration analysis before and after welding to detect 
reduced damping due to the welding process. 
B. DESIGN OF VISCOELASTIC AND CONSTRAINING LAYER 

THICKNESSES 

To approximate system loss factor, the plate structure 


was represented by a three layer beam system and subjected 


to the Bernoulli-Euler type analysis described by M.L. Drake 
and developed by Ross, Kerwin, and Ungar [Ref. 8:pp. 1-6; 
Ref. 9]. The viscoelastic material and constraining layers 
were varied incrementally in the Bernoulli-Euler equations 
(3.1) to (3.3) below, using the FORTRAN program listed in 


Appendix A. 





3 = 
Epi Eyl 
- aL yy 2 2 
EI = . a + ——— i —— + E,H,D + E. He (Ho - D) 
H..-D 
2 paae 31 
+ EH, (H3)-D) - am “5 7 Hg Ee EE eae 
=a 
be her OO eee 3.2 
: ; 
3 
E,H, + 5 +g(E)H, +E,H, +EH,) 
H, + H, +2H, +H G 
1 TH, 2 D 
ii: * ales) D ; and g= 2 (Se 
EHHOK 
where 
E = Young's modulus of elasticity. 
G = shear modulus. 
I = second moment of area. 
H = thickness. 
5 Beam Approximation Length, 2 4 (Peam Approximation Width, 2 
= Plate Length Plate Width 
and 


10 


Subscript 1 refers to the base structure. 
Subscript 2 refers to the damping layer. 
Subscript 3 refers to the outer layer. 


No subscript refers to the composite system. 


The maximum system loss factor and the temperature at which 
it occurred were found for each step. The results of 70 
incremental thickness changes were plotted in Figure 3.3. 
The University of Dayton Research Institute, itil 
suggesting damping design procedures, recommends a system 
loss factor of 0.1 for most applications. (Ref. 10] 
Because Equations (3.1)-(3.3) are only an approximation, the 
desired system loss factor was increased to 0.3 in this 
analysis. From this figure and based on a loss factor of 
-3, the values of the viscoelastic damping layer and 
constraining layer were chosen to be .0625 and .25 inches 


respectively. 


C. VISCOELASTIC MATERIAL 

The damping material applied to all three models tested 
was 3M SJ-201 ISD-112 viscoelastic compound with a thickness 
of 60 mil, supplied with suitable adhesive already applied 
to both sides by the manufacturer, 3M, St. Paul, Minnesota. 
During application of the viscoelastic and constraining 
layers, care was taken to allow air to escape from between 


the layers. Entrapped air can result in severe bond defects 


hall 


giving lower natural frequencies and reduced system loss 
factor. [Ref. 21:3 ppees=95 

ISD-112 was chosen for its excellent loss factor and 
shear modulus characteristics in the zero to 100° F 
temperature range for frequencies below 2,000 Hz. This 
trait is exhibited by the broad loss factor curve of Figure 
3.4 and is highly desirable in many applications. Figure 
3.5 shows a University of Dayton nomogram of the 
experimental determined shear modulus and loss’ factor 
properties for ISD-112. To determine shear modulus and loss 
factor values, locate the material temperature on the top of 
Figure 3.5 and the frequency of interest on the right 
margin. The diagonal line of temperature and horizontal 
line of frequency intersect on the vertical line of reduced 
frequency. Reduced frequency is the frequency value of 
interest multiplied by the temperature of interest and 
divided by the reference temperature of Figure 3.5. Reduced 
frequency and reference temperature are used to convert the 
input values of temperature and frequency to the conditions 
used to plot the loss factor and shear modulus curves. 
(Ref. 12:pp. 47-50] The intersection of reduced frequency 
and the plotted curves give the horizontal line of loss 
factor and shear modulus shown in the left margin. While 
nomograms are useful in comparing characteristics when 
selecting a damping material, precise values are difficult 


to read on the log-log scale and the actual equations.of the 
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nomogram curves were used to give accurate values of loss 
factor and shear modulus in the finite element codes. The 
equations for ISD-112 are given below and used in the 
program listed in Appendix A. [Ref. 9:pp. 6-14] 

Log(FR) = log(F) - 12*(T-TZERO) /(525*T-TZERO) (3.4) 
Log(G) = Log(ML) + (2*Log 10(MROM/ML))/(1+FROM/FR)N (3.5) 


Log(ETA) = Log(ETAFRO) + A*B+D*(1-SQRT(1+A7) )/(C/2) (3.6) 


where the ISD-112 viscoelastic material nomogram constants 


are: 
TZERO = 131 
FROM = 400,000 
MROM = 2059.5 
N = .27 
ML = 28.427 
ETAFRO = 1.25 
SL = 0.27 
SH =—=0:535 
FROL = 150,000 
Cc = PRS, 
and: 
FR = reduced frequency 
F = frequency 
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At = temperature 


G = shear modulus 

ETA = loss factor 

A = (Log(FR) - Log(FROL))/C 
B = SL+SH 

D = SL-SH 


D. NASTRAN FINITE ELEMENT MODEL 

Computer analysis of the sandwich structure was 
performed using the MSC/NASTRAN finite element program. As 
shown in Figures 3.6-3.8 the base plate and constraining 
layer were constructed of QUAD4 plate elements with 3/16 and 
1/8 inch offsets, respectively. The more complicated HEXA8 
solid elements used to simulate the sandwiched viscoelastic 
layer were required to accurately represent the strain 
energy due to shearing of the viscoelastic core and to 
account for the coupling between stretching and bending 
deformations. [Ref. 73:p. 87 Ref. 13¢See.. 2.22-458) The 
Support bulk data card was used to give free-free end 
conditions to the plates. After constructing the finite 
element models, two different analyses were performed to 
predict frequency response amplitude. 

1. Modal Strain Energy (MSE) Method 

Rearranging Equations (2.2) and (2.3) and solving 


for the response amplitude gives the following equation: 


14 


( N oR (wf) *} .op :L)/u," | 2 
tO). = a ee ee 
}¢ 


{l - (w/w)? eat : AS I 
vl 2 ; 
. (5 OO) 5 O/4y 3)? Op, 19R al) \ 2,172 
4 240 My ; 
i ee ing) 5 Wy) 


Vz 
ay (w/a, 5) 





}2 
where: 
X(w) = amplitude of response as a function of 
frequency 
%,,4 = ith natural radian frequency 


°m,i = mode shape value for the excitation grid 
point at the ith natural frequency 


oR, 4 = mode shape value for the response grid 
point at the ith natural frequency 


ns = iterated viscoelastic material loss factor 


(Vy/Vp) 4 = fraction of elastic strain in the visco- 
elastic material at each model 


L = magnitude of the load applied at the 
excitation grid point. 
The MSE method utilizes the MSC/NASTRAN method to find the 
constant values on the right side of Equation (3.7). The 
response amplitude is then calculated external to the finite 
element code. The mode shapes, natural frequencies and 
fractions of elastic strain energy in Equation (3.7) were 
determined using a single MSC/NASTRAN normal mode analysis 
for each plate. The damping material properties input to 
the MSC/NASTRAN normal analysis must be constant and were 
determined using the nomogram Equations (3.4)-(3.6) for a 
temperature of 60° F and the mid-range frequency of 500 Hz. 


The output of the MSC/NASTRAN normal analysis gave the 
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natural frequency, mode shapes, and the fractions of strain 
energy in the damping material at each associated natural 
frequency. The simple plate MSC/NASTRAN normal analysis 
program is listed in Appendix B. 

The loss factor of Equation (3.7) is sensitive to 
the natural frequency and several MSC/NASTRAN mode analyses 
must be applied to compensate for errors in the derived 
natural frequencies. Johnson and Kienholz describe an 
iterative method to approximate this loss factor frequency 
sensitivity. (Ref. 7:pp. 27-31] The iterative method was 
applied by first plotting the damping materials shear 
modulus as a function of frequency using the nomogram 
Equations (3.4)-(3.6). The damping material shear modulus 
of each plate was then varied through several MSC/NASTRAN 
modal analyses and the resulting change in each natural 
frequency was plotted on the same curve. The intersection 
of the curves gives an approximation of the natural 
frequency and shear modulus for each mode. Figure 3.9 
illustrates the iterate technique for two of the simple 
plate modes. The loss factor value was then found from this 
approximate shear modulus using the nomogram Equations 
(3.4)-(3.6). When this iterative method was first applied 
by Johnson and Kienholz, they noted the results differed 
from those observed in frequency response methods. Modes 
below the mid-range frequency show too much damping and 


those above show too little. (Ref. 7:p. 12]. Johnson and 
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Keinholz suggested multiplying the loss factor obtained in 
the described iterative method by (CC/enare mn, where G; is 
the viscoelastic material shear modulus at the ith natural 
frequency and Greer is the core shear modulus at the mid- 
range frequency. Once the matrices of natural frequency, 
mode shape, strain energy fraction and loss factor were 
calculated for each plate, Equation (3.7) was solved and the 
frequency response amplitude was plotted as a function of 
frequency. 
2. Frequency Response Method 

Use of the NMSC/NASTRAN model frequency response 
analysis to directly calculate response amplitude within the 
finite element program requires the viscoelastic material 
properties to be input as constants throughout the frequency 
range analyzed. This restriction limits the modal frequency 
response analysis to frequency spans over which the visco- 
elastic materials shear modulus and loss factor have little 
change. Inspection of the nomogram in Figure 3.5 indicated 
a frequency span of no greater than 100 Hz would give little 
change in the viscoelastic material properties. In order to 
give the best accuracy at the peak response values, the 
natural frequencies of each plate were estimated by 
conducting a MSC/NASTRAN modal analysis using the mid-range 
viscoelastic material properties. The viscoelastic material 
shear modulus and loss factor were then obtained for the 


plate structures at each of these natural frequencies. 
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Modal frequency response analysis was then performed over 
100 Hz spans centered at each natural frequency. The output 
magnitude versus frequency curves for each resonant 
frequency were compiled and plotted as composite curves. 
The modal frequency response analysis in MSC/NASTRAN for the 


simple plate is provided in Appendix B. 


E. VIBRATION GENERATOR TESTING ARRANGEMENT 

Experimentally determined vibration response analysis 
utilized a Wilcoxon research F3 vibration generator as the 
input force device. The F3 vibration generator was attached 
to the base plate on all applications. An HP 3562A Dynamic 
Signal Analyzer provided the swept sine driving signal which 
was amplified by a Wilcoxon research power amplifier (Model 
PA7C). <A force transducer was located before the attachment 
port at the base of the vibration generator to measure the 
input excitation force. The force signal was amplified by a 
charge amplifier and fed as the input signal to the HP 3562A 
Dynamic Signal Analyzer. An Endevco accelerometer was 
placed at a symmetrically opposite point on each plate and 
coupled with a signal conditioner to provide the output 
Signal to the HP 3562A Dynamic Signai Analyzer. Figures 
3.6-3.8 show the locations of the vibration generator and 
accelerometer. The HP 3562A Dynamic Signal Analyzer 
converted the signals to digitized data and computed the 
frequency response, coherence and power spectra. 


Descriptive points were taken from the HP 3562A Dynamic 
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Signal Analyzer display of frequency response and plotted on 


the same graph as the finite element predictions. 


io 





Figure 3.1 Dimensions of Simple Plate with Applied 
Viscoelastic and Constraining Layers 
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Figure 3.2 Dimensions of Milled Plate with Applied 
Viscoelastic and Constraining Layers 


‘a 





szsAeq bututezzsuo\D 
pue Hbutdweq OTASeTecoOsTA FO ssessouyoTyL HutAzen 
U3aTM SReTd oTduts syR FO AOROeT SSOT wa3asAS Jo ajzewtasgA ¢°¢€ eANbTY” 


| AUN IVa d Na 


90°0 


oro 


or de tee  SSANNITHL 3 aoe ae ae aa mae ee ar ee fess 
TOS USAW TO ONTMNVG OF ob bg 
: t- JILSV]IOISTA 


CHBewme et Ftee se M8 oe 


920 3g20 et'0 FTO 
YOLOIVY SSOT WALSAS WAWIXVAN 


0€£°0 


i 3 ssanxgtng 198" 
pesesbeened -MAAWED cbend gt 
‘ONENIVULSNOD : 


Teen 


er he ee eee Bee te ts o 


B8e°0 3=—s« EO 


ev'0 





22 


sinjersedueay pue Aosuenbezq BHbutAzeA YRATM ZTT-dSI jo uoTzeTszep 


; (ZH) AONANOTYA 
oooT OoT 


zojoeI Sssot 


0S 





v 


YOLOVA SSOT 


¢ aanbty 


25 





| 
“/ 
| 


: 
f 
/ 
\ 





TI,» 


AYN 0M UN LA 


Aw i 


ity of Dayton’ Nomogram for 3M ISD-112 Damping Materia 


10 °l 49 


12 


10 


10° 


il 





1/4" Aluminum Constraining 
Layer QUAD4 Elements 


1/16"ISD-112 Viscoelastic 
i Ee 3 
| | Material HEXA Elements 


| yo 


| | 3/8" Aluminum Base Plate 
2.33 IN TYPICAL. QUAD4 Elements 


Figure 3.6 MSC/NASTRAN Finite Element 
Model for Simple Plate 
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Figure 3.7 MSC/NASTRAN Finite Element Model for 
Milled Plate with Free Constraining Layer 
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Figure 3.8 MSC/NASTRAN Finite Element Model for Milled 
Plate with Welded Constraining Layer 
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LV. ea eoULES 


A series of response amplitude versus frequency plots 
were constructed for each plate. The first set of figures 
(Figures. 4.1-4.3) plots the experimental data for each 
plate before and after the application of the viscoelastic 
damping and constraining layers. After completing the two 
prediction methods, the results were plotted (Figures 4.4- 


4.6) with the experimental analysis for the damped plates. 


A. SIMPLE PLATE 

The experimental frequency response plot of the simple 
plate (Figure 4.1) shows a large increase in damping when 
the viscoelastic and constraining layers are applied. All 
natural frequencies shift lower in frequency and in some 
cases amplitude peaks are combined into a single peak. The 
system loss factor for the first mode increased from 2.2% to 
34% and the third mode increased from 1.6% to 14% after the 
application of the viscoelastic material and constraining 
layer. 

Comparison of the two simple plate prediction methods 
with each other and the experimental results show several 
differences (Figure 4.4). The MSE curve resonant 
frequencies differ from the NASTRAN composite curve, with 
the MSE resonant frequencies being slightly higher than the 


NASTRAN composite curve at frequencies below 500 Hz and 


Zo 


slightly lower at frequencies above 500 Hz. The resonant 
frequencies of the experimental data are lower than either 
prediction method throughout the frequency range inspected. 
The experimental data response peak amplitudes are higher 
than either of the prediction curves. The response peak 
amplitudes of the MSE method have approximately the same 
Value as the NASTRAN composite curve with the poorest 
matching above 600 Hz and the best matching at the first 
peak. The experimental data show a loss factor of 34% for 
the first mode while both prediction methods indicate a 45% 
loss factor. The third mode has a 14% loss factor for the 
experimental data and a 20% loss factor for the prediction 


methods. 


B. MILLED PLATE, FREE CONSTRAINING LAYER 

The experimental frequency response plot shows a large 
decrease in peak response amplitudes after the viscoelastic 
damping treatment (Figure 4.2). In several cases, multiple 
amplitude peaks are again combined into a single peak. The 
system loss factor in the 300 Hz range was 1% prior to 
damping. Application of the damping material increases the 
system loss factor to 4% in this range. 

The prediction methods were plotted for the finite 
element model in Figure 4.5 with the two prediction methods 
showing good correlation. At frequencies below 500 Hz both 
prediction methods show similar peak response amplitudes but 


at differing frequencies with the NASTRAN composite method 
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having consistently lower peak frequencies. Above 500 Hz 
both prediction methods again show similar amplitudes but 
the MSE method showed lower peak frequencies. The 
prediction curves are poorly matched above 900 Hz. The 
experimental data show similar resonant frequencies but the 
frequency peaks are sharper than either prediction method 
over the observed frequency range. The experimental system 
loss factor at the 300 Hz peak is 4% while both prediction 


methods indicate a 6% loss factor. 


C. MILLED PLATE WITH WELDED CONSTRAINING LAYER 

After the experimental analysis of the unwelded milled 
plate, the plate constraining layer was welded to the base 
plate while maintaining the majority of the plate surfaces 
below 200° F during the welding process. After welding, the 
aluminum constraining layer and base plate were visibly 
warped and the plate resonated loudly when struck. The 
experimental analysis of the welded plate confirmed a 
Significant decrease in damping (Figure 4.3). The 
experimental loss factors of the welded plate at all natural 
frequencies were indistinguishable from the experimental 
loss factors for the milled plate prior to the application 
of the viscoelastic layer. AFter the weld was cut, the 
plate once again exhibited high damping (Figure 4.7). The 
plot of the two finite element predictions (Figure 4.6) 
yielded curves matched for frequencies below 700 Hz and 


curves beginning to disagree above 700 Hz. 
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V. CONCLUSIONS 


Figure 4.1 shows the simple plate vibration response is 
reduced significantly after the 3M ISD-112 viscoelastic 
material and constraining layer are applied. Comparison of 
the experimental data and the two prediction methods in 
Figure 4.4, however, show that both the NASTRAN composite 
curve and the MSE method predict more damping than has 
occurred. Vibration energy absorbed by the viscoelastic 
layer is converted to heat and possibly increased the 
damping material temperature. No attempt to hold the plate 
temperature constant during the experimental analysis was 
made. Figure 3.9 shows that increasing the temperature of 
the viscoelastic material will decrease its loss factor. 
Another contributing factor to the lower actual damping may 
be the variation in the behavior of the viscoelastic 
material. 

In examining the MSE method (Equation (3.7)), it can be 
seen that the ooff-resonant frequency terms have _= an 
increasing contribution to the predicted response amplitude 
as the frequency increases. 

The natural frequencies of the MSE method are found 
uSing constant mid-range frequency viscoelastic material 
properties. This results in the MSE method amplitude peaks 


of Figure 4.4 shifting downward from the NASTRAN composite 
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curves in the high frequency range. An upward shift of the 
NASTRAN curves occurs in the low frequency range. The two 
prediction methods for the simple plate show more disagree- 
ment than for the other two plate models investigated. The 
MSE summation Equation (3.7) contained 18 modes for each 
plate analysis conducted. However the simple plate is the 
least stiff and thus 14 of these modes occur below 1000 Hz 
(78%). The milled plate had 10 modes below 1000 Hz before 
welding (56%) and 8 after welding (44%). The number of 
terms in the simple plate summation of Equation (3.8) should 
contain more terms to increase the accuracy of the MSE 
method prediction curve. 

Inspection of the experimental analysis in Figure 4.2 
again shows significant damping after application of the 
viscoelastic damping material and constraining layer. 
Figure 4.5 shows that both predicted damping levels were 
greater than the experimental data with the causes of the 
error being the same as for the simple plate. Examination 
of the two prediction curves shows the MSE method gives 
acceptable values of both amplitude and peak response 
frequency throughout the frequency range analyzed. The MSE 
prediction method of the milled plate with ae free 
constraining layer again shows the MSE method response 
frequency shifts described in the simple plate MSE 


prediction. 
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After welding the constraining layer to the base plate, 
the structure was warped and exhibited no damping 
contributable to the viscoelastic material layer (Figure 
4.6). When the weld was cut, the plate warpage decreased 
Significantly and the milled plate once again displayed a 
large loss factor (Figure 4.7). This indicated that the 
viscoelastic layer had separated from the plates as a result 
of the warpage caused by the edge welding process. Because 
of the loss of damping through deformation of the plates 
during welding, only the two prediction methods for the 
milled plate with a welded constraining layer are plotted in 
Figure 4.6. This plate was the stiffest plate analyzed, 
having eight natural frequencies in the frequency range 
studied vice 10 natural frequencies for the milled plate 
prior to welding and 14 for the simple plate. As a result, 
the response amplitude of the MSE prediction method is in 
agreement with the NASTRAN composite prediction method. The 
MSE prediction method again shows a slight upward shift of 
the low frequency range peaks and a downward shift of the 
high range peaks. 

Overall indications of the three plates show actual 
plate behavior is accurately estimated by both prediction 
methods. Errors occurring in the higher natural frequencies 
are due to an insufficient number of terms in the MSE 
summation Equation (3.4). Where inaccuracies occur, the MSE 


prediction method is found to be a good approximation of the 
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NASTRAN prediction method and should be an acceptable 
substitute for the NASTRAN composite prediction method. The 
simplicity of application coupled with the low computer time 
required make the MSE method the most attractive choice for 
engineering analysis which must consider both a wide range 


of environmental conditions and several design configurations. 
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VI. RECOMMENDATIONS 


The following are suggestions for further studies in 


this field: 


1. 


This study investigated prediction accuracies of the 
MSE and composite method at constant temperature. 
Accuracy over a temperature span should be 
investigated. 


The results of the MSE method were heavily dependent 
on the MSC/NASTRAN normal mode analysis conducted 
with viscoelastic properties evaluated at 500 Hz. An 
investigation of the prediction sensitivity to shifts 
in evaluation frequency is recommended. 


The use of the HEXA8 solid finite element to model the 
viscoelastic layer gives results in an MSC/NASTRAN 
program which requires significant computer time to 
solve. Use of the viscoelastic beam element to model 
the viscoelastic layer have been suggested to reduce 
computer time requirements. [Ref. 14:pp. 352-354] 


Failure of the welded plates viscoelastic material was 
due to the high temperature gradient causing warpage, 
yet failure to cool the surfaces of the plate would 
have resulted in destruction of the damping material. 
Investigation of welding techniques is required to 
give practicality to the milled plate design used in 
this paper. 


Application of several viscoelastic material and 
constraining layers have been used on structures in 
order to broaden the frequency range at which damping 
occurs. Applying the prediction methods investigated 
to such multiple damping layers should be researched. 


Viscoelastic material properties may vary from batch 


to batch of the same material. The construction of a 
nomogram uSing samples of the batch applied to the 
structure should greatly increase prediction 
accuracies. 
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APPENDIX A 


NOMOGRAM LOSS FACTOR AND SHEAR MODULUS CALCULATIONS 


As explained in Section III.A, it is impractical to use 
nomogram chart values of loss factor and shear modulus in 
analyzing damping characteristics of systems containing 
Viscoelastic materials. All calculations of shear modulus 
and loss factor in this paper were accomplished through use 
of the FORTRAN program listed in Figures A.1-A.4. Although 
this program is not complicated, use of the equations 
without computer assistance is extremely tedious due to the 
large number of terms. Note that page 4 of the program 
(Figure A.4) also gives the system loss factor based on the 
Bernoulli-Euler equations (3.2)-(3.3). Equations (3.2) and 
(3.3) require Young's Modulus and the shear modulus values 
in order to be solved and thus it is most convenient to 


include these calculations here. 
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APPENDIX B 


MSC/NASTRAN CODE FOR SIMPLE PLATE NORMAL MODE ANALYSIS 


AND _ MSC/NASTRAN CODE FOR SIMPLE PLATE FREQUENCY 
RESPONSE ANALYSIS 

While the MSC/NASTRAN code was used for all the finite 
element plates analyzed, listing each program would serve no 
useful purpose. For completeness, a normal mode analysis 
and a frequency response analysis are listed here for the 
Simple plate structure (Figures B.1-B.10). The MSC/NASTRAN 
programs for the milled plate with free constraining layer 
and milled plate with welded constraining layer vary only in 
the geometric organization of the grid points and finite 


elements. 
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ID DYNAMICS-FREQ MOD 
SOL 3 

TIME 240 

DIAG 849413 

CEND 

TITLESMODAL ANALYSIS S00 HZ PROPERTIES 
ECHO=UNSORT 

SPC=10 

METHOD={] 

ESE=ALL 

DISP=ALL 

OJTPUT( PLOT) 

PLOTTER NASTRAN 
PAPER SIZE 14, BY 10. 
AXES=7,X,Y 

SET 1=sALL 
FIND,SCALE,SET 1,0RIGIN 1 
PLOT 

PLOT WODAL: DEFORMATION 0,1,THRU,15,SET 1,0RIGIN 1,SHAPE 
BEGIN BULK 
GRID0.1e0,0,0,0 

Spx 7)eHeeet(3075) 5% 
=7 
GRI0,2,0,2.33330,0 
Sot (T) op teterX(3.75),= 
=7 
GRID,3,0,4.866/7,0,0 
=ex(7)e=,2,%(3275),5 
=7 

GRI0,4,0,7.20,9 

Spx (7)o=e=5%(3275) 55 
=7 
GRI1D,5,0.69.3533350-9 
See(T)o2eeor(3.75) 5% 
27 
GRID,6,0-11.6667,0,0 
Sex(7) 6 5,2,%(3275) 9 
=7 

GRID,7e001447069 
=ee(7T)ete=e%(32/75),= 
=7 
GRID0,101,0,0,0,.0625 
Sex(7T) gp te=e%*(3275),= 
a7 


Figure B.1 MSC/NASTRAN Normal Mode Analysis of 
the Simple Plate (Page One) 
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GRID, 102,0,2.3333,0,.062e5 
=,*x(7),=,=,%(3./75),= 

=7 
GR1ID,103,0,4.5667,0,.0625 
=,*(7),=,=,%(3.75),5 

=7 

GRID,104,0,7.,07.06eS 
=,%x(7),272,%(3.75),= 

=7 
GRID,105,0,9,333340,.06e5 
s=,x(7),=,2,%(3.75) 42 

=/7 
G2I10,106,0,11.6667,0,.0625 
Sex 7)pse=,%(35.75),= 
=7 
GRI0,107,0,14.,0,.0625 


=ee(7),5556%(3.75),2 

=7 

CQUADG,1,100,1,2,9,8,22°21875 
=pe(b) psp %(7), x07), 807) e077) 0 5e 
=6 
CQIUAD4G,2,100,2,3,104,9542°.,1875 
S,ex(O)pzexC7),xC7T) aC 7) pe T) oot 
=6 
CJUADG,3,100,3,49,11,190,,2°.1875 
=,x CO) pare eC 7) ge C7T) eC 7) ek(7) 6,2 
=6 
CoaUR0G,4,100,4,5,1er,l1lpre-1875 
Spx (6) pte %@(7),eC 7), RCT) oe ( 7) 578 
=6 
CQUAD4,5,100,5,6,13412572°.1875 
=,x(b)p zee C7), R07), eC7T) e077) 66s 
=6 
CQUAD4,6,100,5,7,14413,72°.1875 
Spx (O)o zee) gk(7T) eC 7) oe (7) 55 
=6 
CQUAD4G,101,101-101,102,109,108,,.125 
=px (GO) ep see) geC 7) op k&( 7) pk (C7) oe 
=5 
CQUA04,102,101,102,103,110,109,,.12e5 
Spx (O) pte r(7) eC 7) op ekC7) eo k(7) 5,55 
=6 


Figure B.2 MSC/NASTRAN Normal Mode Analysis of 
the Simple Plate (Page Two) 
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CQUAD4G,103,101,103,104,111,110,,.125 
Sex (Oo aex( 7) px (TF) ,exC7),e(7) 605 

=6 
CQUAD4,104,101,104,105S,112e,111,72125 
Se,x(O)o tek) gx ( 7) ok(7T) ok (7) o ee 

=5 
CQUAD4,105,101,105,106,113,112¢7-2125 
See (O)o see (7) RC 76 eC 7) 6 k( 71) 69 

=6 
CQUAD4G,106,101-106,107,114,113,,e125 
Spt (Ojo Hee) kOT) ROTI G(T) 662 

=§ 

CHEXA,201,200,142,9,8,101-102,401 
Spx(O)oseoX®( 7) gC T) gp kT) oA KC Tp ®%CT) op RCT) A R(1) 
See®(O)e zee ( 7) ge@(T) oe CT) A ® C7) RCT) 0 x07) 6 *%(1) 
Spe (OH) eter (7)pe@CT) A kCT) oe (C7T) A kC To RCT) o 2 C1) 
S,x(O) eter (gk Top kT) o eC) op & C7) ok (7), % C1) 
Sex(O)pzex(T)pe@(7T)ge( Tok (7) ge C7) eo & C7) 6% (1) 
Spt@(O)o tex (7) gp e(T) go eC To RCT) A RCT) oO RCT) A *C1) 
S,x(O)etee(7) RCT) pg RCT) po kC To RCT) A kT) A x C1) 
¢01,109,108 

e{,x7,K%7 

#1 ,%a7,%7 

Kigx7T,x7 

wi,gx7T,x7] 

Ri, x7, x7] 

a1,*7,%7 

a1,a7,nT7 

CHEXA, 202,200,245 3410,9,102,103,+09 

Spt@(O) gp tore ( 7) pk 7T) pk T) GeO A ROT) Go RCT) Go X C1) 
Sprx(O)osok(7T)geC7T) pe (7) po KOT) Gg XC 7) 6 207), % (1) 
Sepx(O) go sor (7) gxl(T)gr(T) op kCT) ok CT) 6 k(7) 9 %*C1) 
Sot (G)op tee (7) pe C7) gp kCT) oe (CT) og kC7T Io C7) 6% (1) 
Sex (O)osor( 7) pg RCT) op X(T) Oo RCT) A hCG BCT) A 2 (1) 


e®(CO) eo Sek ( 7) RCT) Gg RCT) oA RCT) gp kOC7T) 6 (7) 4 x01) 
Sex (Ob) ote eC) ge (C7) op ®( 7) op R(T) op RCT) op &(7) oo &(1) 
409,110,109 

a1,x7,x7 

a1 ,*7,%7 

ei,x],*7 

ai,x7,%7 

Ri ,x7T x7 


Figure B.3 MSC/NASTRAN Normal Mode Analysis of 
the Simple Plate (Page Three) 
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a1 %7,%7 
ai,x],x7 
CHEXA,203,200,3,4,11,10,103,104,+17 
Sex (b)ezee(7) eC 7) ge (7) eC7) ge 7 oe %®(7),2% C1) 
spk (O)o tee 7), eC7) eC 7) op k®( 7) kCT) oe (7) 6% (1) 
=e*(O),=e%(7) R07) R(T) oA RCT) A kK (C7) op kC7T) A 2 OL) 
See (Oo ter (7)ee(7)pxC 7) eC 7) e076 (7) ,%(1) 
=e (CO) ose e(7) RCT) op eC7) op RCT) op R07) op OT) A RON) 
Sek ( Oe Sep e@(7) R(T) R(T) oe (CT) A RCT) op e(7) C1) 
Sok (O)o Soe ( 7) RCT) gp eC 7) go eC TIA RC7) oe (CT) *(1) 
+17,111,110 
k1,*7,%7 
ki,x7,%7 
ai,x/7,x7 
ei,x7,*7 | 
e1,2/7,%7 
ei,x7,x] 
ei,x7,%x7 
CHEXA,204,200,4,5412,114104,105,+2e5 
=pe(b) eae e@(7) eC 7) pe@( To eC 7) pk 7) op &(7) oe C1) 
Spe (O)o ste eC) pe (7) eC T) op eC 7) op eC 7) oe (7),% (1) 
Spe@(O)ep apex (Tg eC 7) op eC To RCT) op RCT) op RCT) A XC1) 
Spe (O)p se e(T),R(C Tp CT) op eC 7) A R(T) op &( 7), * (1) 
See(O)etee(T)gak( 7) pe (7) oe (CT) op k@( 7) ok (7T),%(1) 
Spe (O)paee( Tp kCT) A RCT) oA RITA RCT) A &(7), (1) 
=p kb) op =e R(T) RCT) - R(T) A RCT) Ap KOT) A R(7) 6 *(1) 
25,112,111 

e1L,*]7,%7 
e1,*7,%] 
a1,%x7,%7 
*1,%*7,%7 
a1,x7,%] 
a1,%x7,%/7 
zj1,*x7,*]7 
CHEXA,205,200,5,6,13,12,105,106,+433 
S,ae(b)psere( J) eC 7) ge(7) op e@CT) op 207) oh (7) &(1) 
Spe@(H)o re e( Jee 7) op k@CT) op eC Tp kIT) op eC 7) 5% (1) 
Spe@(O) otek (T)pe( 7) eC 7) 6 ROT) 6 RC 7), e(7),% 01) 
Spe (beter ( 7) pk( 7), R(T) A eC 7g kOCT) A R(T) A RCI) 
See (b)pzex( J) ,e( 7) kCT) oA kC7T) A % C7) 6 207) x C1) 
Sex(B)otex (J )pe@C 7) pkCT) op eC 7), e(7) op &C7) go &O1) 
Spe (bese rkXC2 eC 7) ge (Tee CT k& C7) ok (7) A kC1) 
33,113,112 
el, arias 
e1,27,%]7 

a1L,x/],*] 
ei,xJ,x/] 
k1,%7,*%7 


Figure B.4 MSC/NASTRAN Normal Mode Analysis of 
the Simple Plate (Page Four) 
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£1 ,pe7erl 

ri,gxT,x7 
CHEXA,206,200,6e7414,13,106,107, +441 

=,e (CO) ptr eC) pk T) p(T) op eCT) px C7) ok (7) e201) 


Sek(bH)o tek (7) RCT) ox (CT) ox CT) 6 xC7) oe C7), * (1) 
Sok (BG) oe sor( 7) eO7T) gp eC 7) &C7) A ROT) A kT) A x C1) 
Spe (O)oso rx (7) kT) RCT) A RCT) Ag RCT) A R(T) A C1) 
Sex (G6)eter(7) RCT) gk( To kCT) gp e(T) op k(7) ge (1) 
Spx(H)o te x( 7) px Tg k( Top eC7T) p(T) oe (7) Xx(1) 
Sek (O)p se e®( 7) RCT) gKCT) op eC Top RCT) ok 7T) Ag eC1) 
#41,114,113 

wi,r7T,x7 

ri,x7,x] 

m1,x7,%7 

wi, eT, x7 

RL ,rk7,k7 

a1 ,x7,x7 

rj,x7 x7 

PSHELL,4100,100,.375,100 
PSHELL,101,100,.25,100 

PSOLID,200,200,0 

MATIZ100,1.¢75662334209595 

MAT1,2009 72354-0224 649, .0354651,089 
PARAM,ASINGs 1 

PARAM,AUTOSPC,YES 

SPCi1,10,6,101,THRU,163 

SPC1,10,6-1+/THRU,63 

SUPORT,32,12345 
EIGRelLeAMGIVelee1S00. ge 24e eg tlZ 

#¢Z2Z,MASS 

PARAM,WIMASS, .00259 

ENDDATA 


Figure B.5 MSC/NASTRAN Normal Mode Analysis of 
the Simple Plate (Page Five) 
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ID DYNAMICS,FREQ MOD 
SOL 30 

TIME 240 

CEND 

TITLESFREQUENCY RESPONSE 500 HZ PROPERTIES 
ECHOSUNSORT 

SPC=10 

METHOD] 

SVECTORSALL 
OLOAD=2000 
FREQUENCY=2100 

Se) f£112s24 
DISPLACEMENT (PHASE) =111 
BEGIN BULK 
GRID,1,0,0,0,0 
=,x(7),=,=,%(3.75),= 
=/ 
GRID,2,042.3333,0,0 
=,x(7),2=,5,%(3.75),= 
=/7 

GRID, 3,0,4.6667,0,0 
=,2(07),5,5,%(3.75),2 
=/7 

GRID,4,0,7.70,0 
=,x(7),=,=,%(3./75),2 
=/7 
GR1ID,5,0,9.333370,0 
=,%(7),5,5,%(3.75),= 
=/7 
GR1ID,6,0,11.6667,0,0 
=,x(7),2,=,%(3./75),= 
=7 

GRID,7,0,14.,0,0 
=,%*(7),=,2,%(3.75),5 
=7 
GRID,101,0,0,0,.0625 
=,x(7) 9 =¢=,%(3.75),= 
=7 
GRID,102,0,2.,333370,.0625 
=,*(7),2,=,%(3.,.75),= 
=7 
G210,103,0,%,6667,0,,0625 
=,*(7),=,5,%(3.75),5 
=7 


Figure B.6 MSC/NASTRAN Freguency Response Analysis 


of the Simple Plate (Page One) 
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GamDEaOde Ole 
=,*(7),=,=,% 
=7 
GRI0,105,0,9,333320,.0625 
=,a(7) eter %*(3e75),= 

=/7 

GRID,106,0,11.6667,0,.0625 
=,a(7)525=,%(3.75) 5 

=7 

GRID,107,0,14.,0,.062e5 
Spe(T) pee tet (375) 52 

=7 

CQUAD4,1,100,1,2,9,8,57.-1875 
S,x(b6)paek@(T)gk(T)pe®C7) op & C7) a= 
=6 

CQUAD4,2,100,2, 3,10,96,27.-1875 


Spx (O)etexC7) eC 7) peC 7) oe ( 7) eee 

=6 

CQUAD4,3,100,3,4,11510¢52.1875 

Spx (6) ep zee(7), (7) pg k( 7) pk (7) oe 

=6 

CQIVUADG,4,100,4,S,levlilgr-.1875 
=,R(O)etex( 7) ,x(7),x(T) ox 7) ye 

=6 

CQUAD4,5,100,5,6,13,12,,2°.1875 
=pe(O)e=ee@(7T)gxC 7) pe C 7) pex( 7) eee 

=6 

CQUADG,6,100,5,7,14,13,5°21875 

Spe (OO) op =eX( 7) eC 7) kT) e( 7) 9 o= 

=6 
CQIVUADG,101,101,101,102,109,108,,.125 : 
Zex(6)esex(7),2(7),%(7)px(7)5,2 

=6 
CQUAD4,102,101,102,103,110,109,,.125 
ZSe,x(6)p zee (7 yx( 7) gx 7) pe(7) oe 

=6 
CQUAD4G,103,101,103,104,111,110,,.125 
=SpR®(O)e=ex(7) R(T), R(T) ee (7) 492 

=6 
CQUAD48,104,101,104,105,1l1e,111,,2125 
Spex(6),=ee(T),aCT) op C7) op e(7) eG o= 

=6 
CQUAD4G,105,101,105,106,113,112,.4.12e5 
Spx (HO) gz=ex( 7) gxe( 7) pe@C 7) ex( 7) ge 

=6 


0,.0625 


e828 
(3.75),5 





Figure B.7 MSC/NASTRAN Frequency Response Analysis 
of the Simple Plate (Page Two) 
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CQUAD4G,106,101,106,107,114,113,776125 
=,*(6),2,%(7),%(7),%(7),%(7T) ore 

=6 

Cina geo, ce 00,1,259°767101710e,+0)1 
=,*(O)o=ee(7) eC 7) A RCT) op RCT) op RCT) R(T), %(1) 
Sep*(6) eae eC 7) eC7) eC 7) gp e( 7) RCT) op R(T), *C1) 
=e2*(6),=7%(7),%(7), C7) op e(7) op &C7) op *®(7), (1) 
Sp*(6) p=, %07) e807) p R(T) op RCT) op eC 7) op e(7),*(1) 
=ee(O) pe eC 7) ,R(7) op %( 7) p(T) op %(7) oe %(7) 6 %(1) 
So h®(6) op So k(7) R(T) p R(T) Oo R(T) op RCT) op R(T) eC1) 
Zp%(6),5,%07) 6,807) RCT) op aC 7) 9 eC 7) RCT) op RCL) 
¢01,109,108 

KL ,x7,%7 

R1,%7,%7 

ei,x7,%7 

ae ea a | 

aL ,x7,x7 

wa, 7, «7 

e1,*7,%7 

CHEXA,202,200,2,3,10,9,102,103,#09 

Spe (6), =e%(7) 5807), %07) 6 %07) op (7) oR (7) (1) 
Zp,x(6),5,%(7),k07) pe ( 7) op eC 7) op (7) op RCT) &(1) 
Sp%(6) p56 %(7), R07), C7) op k(C7) gp RCT) ok (7) op a (1) 
Spx(O)pSe%x(7),kC7) pe 7) op eT) op R(T) 6 &(7), (1) 
Sere(6)pa,%8(7),%(7) gp k(7T) AKC 7)» eC 7) op e(7) 4 %(1) 
Sep2(6) parC 7), 807) p kCT) op eROCT) op RCT) 6 ROT) A (1) 
Sep%(6) p Sek 7), KC7) RCT) op RCT) op RC7) op eT) A C1) 
409,110,109 

e1,*7,%7 

eL,*x7,%7 

ei,x7,x/7 


e1,%x7,%7 

ai,x7,x7 

kR1L,27,%x7 

1, %7,%7 
CHEXA,203,200,3,4,11,10,103,104,417 

=e hCG) er Se e(7) eC 7) p(T) eR(7) 4 &(7) (7), RCL) 
Se ® (GO), =e%(7), 807), %(7) p07) 6 807) 6 407) #1) 
=o (6) 54407), 807), %(7),e (7), %(7),e(7) ph (1) 
P*#CO) oar e (7) C7) R(T) R(T) op R(T) op R(T), #1) 
o®(CG) esp %07), 807), %(7) (7) p(T) eo (7), *(1) 
e®CO)e=e%(7),R(7),e (7) o (7), %(7) op &(7) pe(1) 
So ® CO) esp #(7) 807), %(7) (7) pe(7),%(7),8 (1) 
#17,311,110 

ei ,*x7,%7 

kl,x7,x7 

e1,%7,%7 


Figure B.8 MSC/NASTRAN Frequency Response Analysis 
of the Simple Plate (Page Three). 
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Ri,«?,xi 
e1,%x7,x] 
a1,*/7,x7 
a1,*7,x7 
CHEXA,204,200,4,54125,1154104,105,+4+25 
Sex(H)o=zore(T)ge( 7) RCT) op RCT) RCT) oA &(7T) 5 x (1) 
Spx (O) oto x ( 7) gkCT) gk Top X(T) op k(C 7) OA R(T) A *C1) 
Spxe(b)o soe (7) eC T) gp RCT) po &R(T Ap k®CT) A RCT) A k(1) 
Sok(H) ote R(T) Gg kKCT) pA RCT) Go RCT op RCT) o hC7T) A k(1) 
Sot(S)otek(7) p RCT) pkCT) eo e(T) Go hET)oe(7),*%(1) 
SokK(H) ep SeeC7) ROTI gp RCT 6 kCTI ORC TI G6 R(T) oO RCI) 
Spe(bO)ptek®(7T) R(T) RCT) po eOT) op RCT) pk CT) G8 (1) 
42Se,ileslii 
r1,x],%*7 
w1,x7T,x7 
el,x],n] 
ki,x],*7 
e1,*7,%*7 
e1i,x],x] 
S14] , Mt 
CHEXA,205,200,576713,127105,106,+433 
SpRe(O) ete k@CT)pkCT) op e@CT) ok CT eo kCT) oA RCT) A C1) 
Set (O)otore (CT) pkCTI ok CF ok (TI ot CT oe (7) 7% (1) 
Spx(b)egster( 7) ,xC7) gp eC Tok (CT) Ge kCT) ge C7T) A % C1) 
Sprx(H)o tex 7) gp kC7T) pk Tok (CT) A XC T) op kRC7) A *C1) 
SpR(H)o te x(T)gkCT) gp k(T) ok (C To RCT) go kC7) HCI) 
See (CO) otek (7) KC T) gk C To kC TD) pk C7) ok C7) ge C1) 
=,x(H)eptoX(T)pkC 7) px (Tok (Tok CT) ok CT) ok (1) 
#33,113,112 
e1,%7,%7 
a1 ,e7,%7 
ai, %7,%7 
mi,xr7 a7 
e1,x7,xT 
ei ,e7T, 7] 
ei ,x7T,x? 
CHEXA,206,200,6677,14,137106,107,+4+41 
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Figure B.9 MSC/NASTRAN Frequency Response Analysis 
of the Simple Plate (Page Four) 
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Figure B.10 MSC/NASTRAN Frequency Response Analysis 
of the Simple Plate (Page Five) 
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